Improved limit on the permanent electric dipole moment of ^^^Hg 
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We report the results of a new experimental search for a permanent electric dipole moment of 
^''^Hg utilizing a stack of four vapor cells. We find d(i^''Hg) = (0.49±1.29stat±0.76,a^t) x lO"^'' ecm, 
and interpret this as a new upper bound, |d(^^^Hg)| < 3.1x10"^^ e cm (95% C.L.). This result 
improves our previous ^^^Hg limit by a factor of 7, and can be used to set new constraints on CP 
violation in physics beyond the standard model. 

PACS numbers: 11.30.Er,32.10.Dk,32.80.Xx,24.80.+y 



The existence of a finite permanent electric dipole mo- 
ment (EDM) of a particle or atom would violate time 
reversal symmetry (T), and would also imply violation 
of the combined charge conjugation and parity symme- 
try [CP] through the CPT theorem [1 S [3|. EDMs 
are suppressed in the standard model of particle physics 
(SM) , lying many orders of magnitude below current ex- 
perimental sensitivity. However, it is thought that addi- 
tional sources of CP violation are needed to account for 
baryogenesis [1, and many theories beyond the SM, 
such as supersymmetry Q, naturally predict EDMs 
within experimental reach. 

Experimental searches for EDMs have so far yielded 
null results. The most precise and significant limits have 
been set on the EDM of the neutron 0] , the electron Q , 
and the ^^^Hg atom [l3|, leading to tight constraints on 
supersymmetric extensions of the SM [7| . Here we report 
the first result of a new mercury experiment, |d(^^^Hg)| < 
3.1 X 10~^^ e cm (95% C.L.), which improves our previous 
limit by a factor of 7 and provides a yet more exacting 
probe of possible new sources of CP violation. 

i99Hg has a ^5*0 electronic ground state and nuclear 
spin 1/2. An EDM of the ground state atom would point 
along the nuclear spin axis and arise mainly from CP 
violation in the nucleus. We measure the nuclear Larmor 
frequency v given by hv = \2fj,B ± 2dE\, where /i and 
d are the ^^^Hg magnetic and electric dipole moments, 
and B and E are the magnitudes of external magnetic 
and electric fields aligned parallel (+) or antiparallel (— ) 
with each other. The signature for d ^ is thus a shift 
in Larmor frequency when E is reversed relative to B. 

As shown in Fig. [1] our new apparatus uses a stack of 
four spin-polarized Hg vapor cells in a common _B-field. 
The middle two cells have oppositely directed i?-fields, re- 
sulting in EDM-sensitive Larmor shifts of opposite sign; 
the outer two cells, enclosed by the high voltage (HV) 
electrodes and thus placed at = 0, are free of EDM 
effects and serve to cancel i3-field gradient noise and pro- 
vide checks for spurious HV-correlated B-field shifts. 

The vapor cells are constructed from high purity fused 
silica and contain isotopically enriched ^^^Hg (92 %) at a 
density of 4 x 10^'^ cm"'^, a paraffin wall coating, and 475 
Torr of CO buffer gas. CO efficiently quenches excited 



state ^^^Hg and thus reduces degradation of the wall 
coating [llj. Spin coherence times T2 are 100 to 200 sec. 
A conductive SnO coating on the cell end-caps provides 
electric field plates separated by 11 mm. The average 
leakage currents across the cells are 0.42 pA at ±10 kV. 

The vapor cells are housed inside a conductive ves- 
sel. Upper and lower feedthroughs connect the elec- 
trodes containing the outer cells to a HV supply. Holes 
in the electrodes provide optical access. All materials are 
free of measurable magnetic impurities. The vessel and 
electrodes are constructed from electrically conductive, 
graphite-filled polyethylene. A gold coated, fused-silica 
plate creates a groundplane between the two middle cells. 
To minimize leakage currents, the vessel is either period- 
ically or continuously flushed with SFg or N2. 

Light for polarizing the ^^^Hg spins and detecting their 



precession is generated with a 254 nm laser system [12 1 
and directed through each cell with the fc- vector perpen- 
dicular to the precession axis along B. During the 30 sec 
pump phase, the light is circularly polarized, tuned to the 
center of the ^^^Hg ^SaiF ^ 1/2) ^Pi{F ^ 1/2) tran- 
sition, and amplitude modulated at the 16 Hz Larmor 
frequency set by the 22 mC main magnetic field, thereby 
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FIG. 1: (Color online) Simplified diagram of the "^Hg EDM 
apparatus showing details of the vapor cell holding vessel and 
middle cell light beams. The topmost cell is shown inside a 
cutaway view of the top electrode, while the bottom electrode 
shows a light access hole for the enclosed cell. The outer cell 
light beams are not shown but travel along the magnetic shield 
axis, perpendicular to the middle cell beams. 
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FIG. 2: (a) Typical single-cell precession signal with an ex- 
panded 0.5 sec segment, (b) and (c) show and Auedm for a 
typical run. In (c) the reduced is 1-2 and the run-averaged 



statistical error is 0.85 nHz after scaling by ^yx' 



building up precessing ^^^Hg spin polarization a by syn- 
chronous optical pumping. During the probe phase, the 
Ught polarization is switched to linear and the frequency 
is tuned midway between the F — 1/2 and 3/2 hyperfine 
lines, a detuning that gives large optical rotation angles 
(proportional to cf ■ k), vanishing circular dichroism, and 
relatively low absorption. The precession of a modulates 
the light polarization angle at the Larmor frequency; the 
angle is measured, for each cell, by a photodiode after a 
linear polarizer. The spin precession is monitored for 
100-200 sec, after which the pump/probe cycle is re- 
peated. The HV is ramped to a new value during the 
pump phase, typically alternating between ±10 kV. 

The Larmor frequencies and their errors are extracted 
by fitting exponentially decaying sine waves to the photo- 
diode signals. Linear combinations of the fitted frequen- 
cies are then constructed. The combination that maxi- 
mally suppresses magnetic gradient noise, and thus gives 
the best EDM sensitivity, is Vc = {^'mt — ^mb) — \{vot^ 
^ob), where the subscripts denote cell positions: OT is 
outer-top, MB is middle-bottom, etc. The factor of ^ re- 
sults from the three times greater separation between the 
two outer cells compared to the middle cells. Combina- 
tions with zero EDM sensitivity are simultaneously used 
to look for systematic effects. The EDM signal Ai^edm 
is obtained from the HV-correlated component of i^c via 
S^oint string analysis that removes linear drifts in time 
[l3| . Figure 2 shows a single-cell precession signal for an 
individual scan, along with Vc and Ai'edm for a typical 
run. Data runs lasted roughly 24 hours and consisted 
of several hundred individual scans. The run-averaged 
statistical error for Avedm is set by the weighted error 
of the mean multiplied by the square root of the reduced 
(for typical runs, was 2 or less). 

Key components were periodically changed or reversed, 
with data taken for nine vapor cells, four electrodes, two 



vessels, multiple vapor cell and electrode orientations, 
and various configurations of the photodiode data acqui- 
sition (DAQ) channels. The vapor cell, electrode, and 
vessel flips used nominally identical components. Com- 
ponent changes were made between groups of 10-20 runs 
termed sequences; in parallel, the paraffin inside each 
cell was remelted and the outer surfaces of the cells were 
cleaned. Each sequence comprised a roughly equal num- 
ber of dipole HV runs (H 1 — HV sequence) for the 

two main _B-field directions; one or more runs with a 
quadrupole HV sequence (0-1-0— 0-f), sensitive to ef- 
fects; and several tilted field (TF) dipole HV runs with 
the main _B-field tipped by ±10° along the middle and 
outer light beam A:- vectors, sensitive to v x E motional 
B-fields. A limited number of dipole HV runs were taken 
at 7 kV and 5 kV. The TF runs and two sequences at high 
light intensity (used to set limits on intensity-dependent 
shifts) were excluded from the final value for Avedm due 
to susceptibility to additional systematic errors. 

An unknown, HV-correlated, EDM-mimicking offset 
was added to the fitted values of vmt and vmb- This 
fixed blind offset masked the measured EDM and was re- 
vealed only after the data collection, data cuts, and error 
analysis were complete. 

A potential correlation (> 90% probability) was found 
between Avedm and the number of micro-sparks per 
scan. Nearly every spark, logged as short duration, > 
100 pA spikes in the continuously monitored leakage cur- 
rents, occurred in 5 sequences when the vessel was peri- 
odically flushed with N2. Two approaches to removing 
this apparent correlation were tested: one cut individ- 
ual scans where sparks occurred, the other cut entire se- 
quences containing scans with sparks The resulting cen- 
tral values for Avedu in the two approaches agreed to 
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FIG. 3: ^^^Hg EDM versus sequence number. Open (closed) 
symbols denote +_B {—B). Squares (circles) denote vessel 1 
(vessel 2). The lower-left inset is a histogram of d(^^^Hg) 
for the 166 separate runs. The lower-right inset shows the 
dataset-wide -\-B and —B values along with the final ci(^^^Hg). 
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within 1.8x 10"^" e cm. The change to I^vedm with and 
without spark cuts was < 4.2 x 10~^° e cm. We adopted 
the more conservative sequence-ehmination approach. 

Figure [3] shows (i(^^^Hg) for the 15 sequences that 
passed the selection criteria, along with a histogram of 
(i(^^^Hg) for the corresponding 166 runs. The sequence 
values are divided into one point for each i?-field direc- 
tion; each point is the weighted average of the relevant 
runs within the sequence. In each case, the +_B and —B 
data are in good agreement. As shown in the lower right 
inset, the weighted average of all the and —B data 
also agree within 1-cr. Systematics that change sign (rela- 
tive to the EDM signal) when B is reversed would appear 
in the difference, but cancel in the average of the -fS 
and —B results. Although the data is apparently free 
of such problems, we determined sequence-level values 
from straight averages of the ^B and —B data. The cen- 
tral value for the entire dataset was then obtained from 
the weighted average of the sequence values: (i(^^^Hg) = 
(0.49 ± Xfl^stat) X 10~^^ e cm. The statistical error corre- 
sponds to a frequency difference between the two middle 
cells of 0.1 nHz, a 4x improvement on Ref. (lo| . 

With individual run errors set as discussed above, the 
reduced for the 166 runs is 0.65. When the individual 
scans are grouped into 3 hour segments and the same 
procedure is applied, ^ \. One potential source of 
this behavior is low frequency drift that averages faster 
than white noise, due for example to beam pointing drift 
that is tied to the 3 hour time scale for resets of the 
piezo-actuated laser cavity mirrors. 

Figure |4] summarizes several checks for systematic ef- 
fects. We did not find statistically significant correlations 
between ^vedm and the vapor cells or electrodes (or 
their orientation inside the vessel), the DAQ channel or- 
dering, or the vessels. Values for d(^^^Hg) extracted from 
the TF runs, the quadrupole HV runs, the high intensity 
runs, and the two-cell difference, MD = {vmt — i^mb), 
agree with the Fig. [3] final value at the l-cr level. 

Table U summarizes the systematic errors. Three con- 
tributions dominate. The spark analysis error is the 
difference in the final Avedm value with and without 
spark cuts. The parameter correlations error is obtained 
by multiplying the HV correlation for key experimen- 
tal parameters by the correlation of each parameter with 
Avedm, and then summing the products in quadrature. 
Specific parameters are: the vapor cell spin amplitudes, 
lifetimes, relative phases, and light transmission; the 
laser power, frequency, drive current, and piezo control 
voltages; an external 3-axis fluxgate magnetometer; and 
the -B-field coil currents (main coil and 3 gradient coils). 
No statistically significant correlations were found. The 
leakage current induced error cannot be obtained from 
the (unresolved) correlation slope oi Avedm versus leak- 
age current due to the limited range over which the cur- 
rents varied. A conservative estimate is instead obtained 
from the worst case scenario of current flow along heli- 
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FIG. 4: Dependence of d(^^^Hg) on the vapor cells, electrodes, 
DAQ channel order, and vessels. The second pane from the 
right shows d(^^''Hg) extracted from the quadrupole HV, high 
intensity, and TF runs and the two-cell difference (MD) . The 
right-most pane shows the final d(^^^Hg) described in the text. 

cal paths within or on the vapor cell walls. The vapor 
cell geometry, with two opposed seal-off stems, limits a 
helical current path to < 1/2 full turn around the cell. 
The average single-cell leakage current was 0.42 pA. The 
resulting fields in the two middle cells can either add or 
subtract, so we take x 0.42 pA = 0.59 pA as the ef- 
fective current. The measurement used 9 different vapor 
cells (4 dominate the statistical error) whose helical cur- 
rent paths should be uncorrelated; to account for this 
averaging, we divide by 2. Combining the above gives 
a systematic error of 4.53 x 10~^° e cm. The remaining 
Table U entries will be discussed in a longer publication. 

Summing the systematic errors in quadrature leads to 
our final result: 

d(199Hg) = (0.49 ± l.2^stat ± G-lQsyst) X 10^29^ 

which we interpret as an upper limit of |d(^^^Hg)| < 
3.1x10^^^ ecm (95% C.L.). This new limit can be trans- 
lated into upper bounds on more fundamental CP vio- 
lating parameters. Table [Til summarizes these limits. 

CP violation in atomic nuclei is conventionally pa- 
rameterized by the Schiff moment 5, the lowest or- 
der CP violating nuclear moment unscreened by the 
electron cloud. Atomic calculations give (i(^^^Hg) = 
-2.8 X 10-"(S'/efm^) e cm The dominant contri- 
bution to the i^^Hg Schiff moment is expected to come 



TABLE I: Systematic error budget (10"^° e cm). 



Source 


Error 


Source 


Error 


Leakage Currents 


4.53 


Charging Currents 


0.40 


Parameter Correlations 


4.31 


Convection 


0.36 


Spark Analysis 


4.16 


{v X E) B-Fields 


0.18 


Stark Interference 


1.09 


Berry's Phase 


0.18 


E'^ Effects 


0.62 


Quadrature Sum 


7.63 
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TABLE II: Limits on CP violating parameters based on our 
new experimental limit for d(^^^Hg) (95% C.L.) compared to 
limits from the Tl (90% C.L.) [sf, neutron (90% C.L.) 
or TIF (95% C.L.) [13] experiments. Values that improve 
upon (complement) previous limits appear above (below) the 
horizontal line. Relevant theory references for the alternate 
limits are given in the last column. 

Parameter ^^^Hg bound Hg theory Best alternate limit 



rf,(cm) " 


6 xlO" 


-27 


[15j 


n; 


3 XlO" 


26 


[3] 


dp{e cm) 


7.9 xlO' 


-25 


[16] 


TIF: 


6 XlO" 


23 


fl7] 


Cs 


5.2 xlO" 


-8 


[18] 


Tl: 


2.4 XlO" 


7 


[19] 


Cp 


5.1 xlO" 


-7 


[18] 


TIF: 


3 xlO" 
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Ct 


1.5 xlO" 


-9 


[18] 


TIF: 


4.5 xlO" 


7 


m 


Oqcd 


3 xlO" 


-10 




n: 


1 XlO" 


10 


[3] 


dn(e cm) 


5.8 XlO" 


-26 




n: 


2.9 XlO" 


26 


[3] 


de(e cm) 


3 xlO" 


-27 


[21, 22] 


Tl: 


1.6 XlO" 


27 
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n; dq = 
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from CP violating nucleon-nucleon interactions. The 
most recent calculation gives 5'(^^^Hg) = g{0.01g''^^ + 
0.075(1) + 0.02g(2))efm^ [H, where g is the strong 
pion-nucleon-nucleon {nNN) coupling constant, and g 
denotes CP-odd nNN isoscalar (0), isovector (1), and 
isotensor (2) couplings. Note that recent work con- 
siders whether the form of S used to interpret atomic 
EDMs should be modified 



25, \2&i. The isovector cou- 



pling has been calculated in terms of chromo-EDMs of 
the quarks, giving g^^^ = 2(d„ — dd)lO^^ cm~^ [15|. For 
comparison, the neutron EDM is dependent on a dif- 
ferent combination of quark EDMs and chromo-EDMs, 
dn = l.le(0.5d„ -I- dd) + lA{dd - 0.25d„) ^. The quark 
EDMs and chromo-EDMs can be estimated in supersym- 
metric extensions of the SM and in many cases, our result 
for d(i^^Hg) provides the most stringent constraints on 
CP violating phases 0, H, 01- The isoscalar couphng is 
dependent on the CP violating phase in the QCD La- 



grangian, 5^"^ ~ 0.027 120|. 

Although expected to be smaller than the above effects, 
the proton and neutron EDM contributions to S were 
calculated in giving Si^^^Rg) = (l.M„ + 0.2dp) fm^. 
This allows us to set limits, < 5.8 x 10~^^ e cm and 
I dp I < 7.9 X 10~^^ e cm, where in the case of t he p roton 
we include a 30% theoretical uncertainty as in nj. This 
improves the best upper limit on dp by a factor of 7. 

The i^^Hg EDM can also have possible contributions 
from CP violating semileptonic interactions between the 
atomic electrons and nucleons, typically parameterized 
in terms of the constants Cs, Cp, and Ct 0, New 
bounds on these parameters are shown in Table [Til We 
also present a limit on the electron EDM derived from 
the hyperfine structure coupling between the nuclear and 



In summary, we have performed a new search for the 
i99Hg EDM. We improve the previous limit by a factor 
of 7, placing new bounds on hadronic and semileptonic 
CP violation. We are currently upgrading the apparatus 
and expect further improvements in the EDM sensitivity. 
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David Meyer, Bob Morley, and Kristian Waedegaard. 



electronic spins, d(i^^Hg) ~ 0.012 4 [2^. This hmit, 
however, is relatively uncertain since an alternate esti- 
mate is of opposite sign, di^^^Rg) ~ -0.014 4 U^^- 
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